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BACKGROUND 


Magnetic  field  gradient  NMR  is  well  adapted  to  measuring  diffusion 
in  colloidal  systems,  including  biological  cells,  because  the  experi- 
mental measurement  times  by  this  method  are  such  that  the  distances 
traveled  by  the  molecules  are  of  the  same  order  as  the  dimensions 
of  the  inhomogeneities  of  the  system.  The  result  is  that  the  apparent 
diffusion  coefficients  are  dependent  on  the  diffusion  time.  By  vary- 
ing this  latter  parameter  the  dimensions  of  the  inhomogeneities  as 
well  as  the  local  diffusion  coefficients  within  them  can  be  obtained. 

The  maximum  information  is  obtained  by  the  use  of  the  widest 
possible  range  of  diffusion  times.  In  previously  reported  work,  a 
variety  of  recently  developed  techniques  have  been  used  to  get  a 
much  wider  range  of  diffusion  times  than  had  been  employed  in  earlier 
studies  of  diffusion  in  biological  materials. 

The  theory  of  the  use  of  pulsed-field-gradient  NMR  to  measure 
intracellular  self  diffusion  has  been  outlined  in  Annual  Report  No. 

I1  of  this  project,  along  with  a description  of  the  equipment  assembled 
to  perform  the  measurements.  A preliminary  description  of  the  results 
of  a series  of  measurements  on  various  frog  muscles  was  also  given 
there. 

The  experimental  results  were  presented  in  full  in  Annual  Report 
No.  2 of  this  project,  including  diffusion  and  spin  relaxation 
measurements  on  frog  muscles,  red  cells,  yeast  and  E.  coli.  They 
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are  presently  being  written  for  publication  ip  the  open  literature. 

Much  of  this  work  was  presented  at  the  21it  Annual  Meeting  of  the 
Biophysical  Society,  February  15-18,  1977  in  New  Orleans. 

Since  then,  a complete  theoretical  derivation  of  the  relation 
between  diffusion  coefficients  and  diffusion  time,  as  measured  by 
an  NMR  pul sed-gradient  experiment,  has  been  performed.  The  results 
constitute  the  present  report.  This  work  has  been  accepted  for 
publication  by  the  Journal  of  Chemical  Physics. 
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INTRODUCTION 


When  diffusional  motion  is  observed  in  an  inhomogeneous  system 
the  degree  of  mobility  of  any  particular  substance  appears  to  depend 
on  the  length  of  the  observation  time.  At  short  enough  times  a wide 
spectrum  of  mobilities  may  be  observed,  corresponding  to  a variation 
in  local  viscosities.  At  long  enough  observation  times  a uniform 
averaged  mobility  is  observed  because  the  substance  has  an  opportunity 
to  repeatedly  sample  all  of  the  different  environments. 

One  particular  type  of  system  which  occurs  frequently  is  a set 
of  barriers  in  an  otherwise  homogeneous  medium.  The  amount  of  sub- 
stance within  the  barrier  at  any  one  time  may  be  quite  small.  In 
that  case,  if  diffusion  is  observed  over  a short  enough  time,  very 
little  of  the  substance  experiences  the  effect  of  the  barriers,  and 
the  observed  motion  is  characteristic  of  the  medium  alone.  As  the 
time  of  observation  is  extended,  more  of  the  substance  is  reflected 
at  barriers;  and  thus,  its  total  displacement  is  less  than  would 
have  been  the  case  without  the  barriers.  We  say  that  the  apparent 
(transient)  diffusion  coefficient  D(tj  decreases,  although  the  total 
displacement  generally  has  increased.  In  the  limit  of  long  obser- 
vation times  the  apparent  diffusion  coefficient  asymptotically 
approaches  a new  value  which  depends  on  the  permeability  and  the 
geometric  arrangement  of  the  barriers  as  well  as  the  viscosity  of 
the  medium. 


5 


Magnetic-field-gradient,  spin-echo  methods  have  been  increasingly 
employed  for  measurements  in  this  type  of  .ystem  when  the  spacing 
of  the  barriers  corresponds  to  colloidal  dimensions  (for  reviews, 
see  references  3 and  4).  This  is  because  the  experimentally 
accessible  range  of  diffusion  times  corresponds  to  distances  of  the 
order  of  colloidal  dimensions,  so  that  diffusion  can  be  observed 
over  a range  of  times  such  that  a considerable  change  in  apparent 
diffusion  coefficient  may  be  observed. 

By  fitting  a set  of  such  diffusion  measurements  to  mathematical 
functions  appropriate  to  the  geometry  of  the  system  under  observation 
it  is  possible  to  extract  the  system  variables,  i.e.,  the  local 
diffusion  coefficients,  and  the  barrier  spacings  and  permeabilities. 

The  appropriate  mathematical  functions  are  usually  quite  com- 
plicated even  for  simple  geometries,  however.  Derivations  are  avail- 
able only  for  homogeneous  media  bounded  by  spherical,^’®  cylindrical,^'^ 
5 8-11 

or  planar  ’ barriers  which  are  impermeable.  The  author  is  aware 
of  no  formula  for  any  geometry  of  periodic  barriers  of  arbitrary 
permeability.  Nevertheless,  such  geometries  are  present  in  many 
interesting  experimental  systems.  The  outer  membrane  of  most  cells 
is  highly  permeable  to  water  and  to  many  other  substances.  The  con- 
tinuous phase  of  an  emulsion  and  open  pore  rocks,  minerals,  and 
plastic  foams  are  also  systems  containing  barriers  which  may  be 
penetrated  or  bypassed. 
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A few  suggestions  have  been  made  of  approximate  methods  xor 
handling  data  from  such  systems,  based  on  a superposition  of  the 
short  time  and  long  time  limits  of  the  diffusion  coefficient,  D(o) 
and  D(<«)  respectively,  according  to  an  equation 

D(t)  « f(t)*  D(o)  + [l-f(t)].  D(»)  (1) 

4 12 

where  f(t)  is  assumed  independent  of  D (<*>).  Derivations  such  as 
in  references  5-11  give  f(t)  for  D(«>)  = 0.  We  may  ask  whether  the 
same  f(t)  cculd  be  used  in  Eq.(l)  when  D(«°)  f 0. 

We  report  here  an  exact  derivation  of  D(jL)  for  one  special  case, 
and  we  show  that  it  cannot  be  approximated  as  in  Eq.(l).  The  chosen 
case  is  one-dimensional  diffusion  through  equally  spaced  barriers  of 
equal  but  arbitrary  permeability,  £.  A finite,  constant  diffusion 
coefficient  D(o)  exists  in  the  intervening  spaces.  In  three 
dimensional  space  the  barriers  would  consist  of  parallel  planes 
normal  to  the  direction  in  which  the  diffusional  observations  are  made, 
see  Figure  1.  Non-absorbing  barriers  parallel  to  the  direction  of 
observation  would  have  no  effect  on  motion  in  this  direction;  therefore, 
the  geometry  chosen  can  also  represent  cubical  barriers.  It  also  may 
be  expected  to  approximate  to  first  order  the  geometry  of  many  practi- 
cal systems  of  regularly  spaced  barriers. 
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The  results  obtained  below  apply  equally  well  to  the  ease  of 
transient  heat  flow  in  a stack  of  planar  slabs  with  finite  contact 
resistance  at  the  interfaces  and  a planar  source.  In  that  case, 
which  we  will  not  discuss  further,  one  should  replace  D(o)  by  « = k /c, 
£ by  h/c,  and  p by  the  temperature,  where  * is  the  thermal  diffusivity, 
J<  is  the  thermal  conductivity,  £ is  the  heat  capacity  per  unit  volume, 
h is  the  "radiation  constant"  at  the  interface,  and  p is  the  spin 
distribution  function. 


DERIVATION 

Our  purpose  is  to  calculate  the  diffusion  coefficient  in  the 

chosen  system  as  it  would  be  measured  by  a spin-echo  NMR  experiment 

where  the  echo  is  attenuated  by  the  application  of  pairs  of  field 

gradient  pulses  of  magnitude  £,  and  of  duration  <5,  short  compared  to 

iheir  separation,  a.  It  has  been  shown  that  in  such  an  experiment 

a represents  the  diffusion  time,  t,  and  that  the  relative  echo 

6 8 

height,  R,  is  given  by  ’ 


R(t)  = 


p(Xq|x.  t)  • cos[Y6g(x-x0)]dxdxo, 


(2) 


where  ()  is  the  conditional  probability  for  the  arrival  at  x,  at  time 
t of  nuclei  which  originated  at  at  t = 0.  R is  defined  as  the 
rftio  of  signal  (echo)  strengths  with  and  without  the  gradient,  and 
is  less  than  unity  if  diffusion  is  taking  place.  Y is  the  nuclear 
gyromagnetic  ratio. 
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or  reguiariy  spacea  pianar  Darners.  The  location  of  the  instantaneous 
source,  Xq,  is  measured  from  the  center  of  the  region  2=0;  whereas  x is 
measured  from  the  center  of  the  region  under  consideration. 


In  a homogeneous  system  the  diffusion  coefficient  is  obtained 
from' ^ 


■>  y '* 

lnR(t)  = -y  Da  <S  t.  (3) 

Assuming  no  prior  knowledge  of  the  inhomogeneity  of  our  system  we 
combine  Eq.(2)  and  Eq . ( 3 ) to  define  a transient  "apparent"  diffusion 
coefficient,  D(t). 

It  is  also  informative  to  define  the  transient  diffusion  coefficient 

in  a different  way,  by  analogy  to  the  Einstein  relation,  as 
2 

(x-j^)  = 2D( t)  • t.  Averaging  over  we  then  have: 

r°°  2 

D(t)  = (1/2 1)^  J"  p (2^1  x dxd*0  (4) 


Both  Eq.  (3)  and  Eq.  (4)  will  be  used,  and  the  results  compared. 

For  a theoretical  derivation  of  either  case  we  need  an  expression 
for  p.  Formulas  from  the  literature,'4’  valid  for  composite 
systems  with  a finite  number  of  layers,  have  been  adapted  to  our 
case  of  interfaces  (barriers)  of  equal  spacing  and  permeability. 
Provided  is  within  the  system,  and  t is  small  enough  so  that  p 
is  negligible  in  the  outer  layers,  the  numerical  results  should  be 
the  same  as  for  an  infinite  number  of  layers.  The  resulting 
distribution  function  is: 


n=i 


E / 

all  i nth 


exp(-Br>D(o)  -t) 


X i„(x)dx 


- in'- 


(f) 
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2 2 2 . 

Here  Xin  are  an  infinite  set  of  eigenfunctions  of  (9  /9x  +6n)  *in  = 0, 
discontinuous  at  the  barriers,  and  of  the  form  = c^coss^  + d^sing^. 
The  layers  are  numbered  by  X (see  Figure  (1)),  and  the  subscript  "o" 
refers  to  the  layer  in  which  2^  occurs.  The  functions  satisfy  the 
orthogonality  relation  V''  f X-  • X.  dx  = 0,  n / m. 

J i th  “"-111  ~ - 

al  1 i 

The  boundary  conditions  are  that  for  eich  interface 

;)X.j/9x  =9Xi.  + )/9x  and  (6) 

£•  (yi+i-X.i  ) = D(o)-3X.  /ax  (7) 

At  the  outer  faces  we  arbitrarily  assume  diffusion  into  a sink, 
so  that  one  or  the  other  of  the  X^  on  the  left  side  of  Eq.(7)  is  zero, 
and  Eq.(6)  does  not  apply.  These  boundary  conditions  form  a set 
of  homogeneous  simultaneous  equations  from  which  sets  of  c,  d,  and  6 
are  obtained. 

For  substituting  Eq.(5)  into  Eqs.(2)  or  (4),  a great  simpli- 
fication of  the  algebra  and  of  the  numerical  computations  was  accom- 
plished by  defining  the  coordinate  so  as  to  let  x represent  the 
distance  from  the  middle  of  the  same  layer,  in  units  of  a,  as  in 
Figure  1.  The  term  U-Xq)  in  Eqs . ( 2)  and  (4)  then  became  (x+1-Xq), 
where  i_  is  the  number  of  layers  from  the  one  containing  to  the  one 
containing  x. 


The  integration  over  x was  carried  out  for  each  layer  separately, 
then  summed  over  all  layers.  If  a large  enough  number  of  layers  is 
chosen  so  that  p ( |_x  * t_)  is  negligible  in  the  outermost  layers  at 
the  longest  times,  t^,  to  be  considered,  then  the  system  has  transla- 
tional symmetry,  and  the  integration  over  need  be  performed  only 
over  one  layer,  for  convenience  the  center  layer  of  an  odd-numbered 
stack. 

The  final  expression  for  the  spin  echo  height  is  then: 

R( t)  = ^ ^ (Z+YHZ-Y)  E (c^.-c.d^)  sinei 

+ W)  ^ E cosqJ_  + (Z-Y)  ^ E d^cosoi  ■]n/V 
exp(-(<nI)  (8) 


where  Z = s inC(il+»n)/2]  f y = sin[(tVBn)/2]  , ©=Y6ga , j = D(o)t/a?,  and 


0+6 


0-6. 


= (1/2)  ^(c.+d^  + (l/Bn)sin(Bn/2)cos(en/2)  S(c^-dJ). 


D(t)  is  then  calculated  from  Eq.(3). 


The  alternate  expression  for  D(t)  as  defined  by  Eg. (4)  is: 


E 


0(t)  = (1/21)  Z^(V%)exP(-6nD* 

n — - — 


(9) 


12 


where  is  as  before,  and 

v <v2y  si"v4sinv2(s,nvs2cos^,]^  % 

-(vv  sinysinvv°sv^i  i 

+<VV  s1n  4,?  £4 
_ _2~ 

-(d0/2bn)[sinbn-bncosbfl]‘  E di  , 

3 „ 

+ (c^/^)sin^[sinb„-^,cosb  ] 2-i  d . i , 

-o  — n — nL  — n -n  — nJ  — v- 

where  t>  = S J2. 

— n n 

NUMERICAL  RESULTS 

Using  Eqs.  (8)  and  (9)  a large  number  of  numerical  evaluations 
were  performed,  mostly  for  a 19-region  system,  although  a smaller 
number  of  computations  were  also  performed  for  1-,  3-,  5-,  7-,  and  41- 
region  systems.  All  of  the  results  presented  below  pertain  to  19 
layers,  unless  otherwise  noted. 
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The  Distribution  Function,  p ( Xq | x_,  t) 

Numerical  calculations  of  p are  not  necessary  to  the  calculation 
of  apparent  diffusion  coefficients  since  these  may  be  obtained 
directly  from  Eq.  (8)  or  Eq.  (9).  Nevertheless,  these  distribution 
functions  are  interesting,  and  their  evaluation  afforded  checks  on 
a part  of  the  procedure. 

Distributions  for  an  instantaneous  source  at  the  center  (x  =0) 

— o 

of  one  region  are  presented  in  Figure  2 for  different  values  of  the 

2 

reduced  permeability,  P = aj)/D(o)  and  of  the  reduced  time,  T = D(o)*t/a  . 
Since  the  functions  are  symmetric  about  x^  The  curves  are  reproduced 
only  on  one  side  of  the  region  containing  the  source.  One  example 
is  given  for  an  assymetrical ly  placed  x^  in  Figure  3. 

The  results  are  in  many  ways  as  expected.  The  functions  decrease 
monotonical ly  on  both  sides  of  a single  maximum,  and  are  discontinuous 
at  the  barriers. 

The  maximum  of  p remains  at  the  center  for  x =0,  or  moves 

— o 

eventually  toward  the  center  if  x^O.  A number  of  numerical  calcu- 
lations of  the  slopes  on  both  sides  of  an  interface  showed  that  Eqs. 

(6)  and  (7)  are  satisfied;  and  calculations  of  p at  several  closely 

? 2 

spaced  T showed  that  6 p / 6_x  = 6p/<sT  is  satisfied  in  the  regions  between 
the  barriers.  Each  of  these  checks  was  done  for  several  values  of 
P and  of  T. 
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£(x,T  ) 


Fig.  2(a).  Distribution  of  randomly  diffusing  particles  from  an 
instantaneous  source  at  the  center  of  the  middle 
layer  of  a 19-layer  system  at  reduced  times  T = 
•1,  1,  10,  as  labeled,  at  reduced  permeability  P 


LAYER 

Fig.  2(b).  Distribution  of  randomly  diffusing  particles  from  an 
instantaneous  source  at  the  center  of  the  middle 
layer  of  a 19-layer  system  at  reduced  times  T = .01, 

.1,  1,  10,  as  labeled,  at  reduced  permeability  P = 2.0. 


p ( X ,T  ) 


Fig.  2(c).  Distribution  of  randomly  diffusing  particles  from  an 
instantaneous  source  at  the  center  of  the  middle 
layer  of  a 19-layer  system  at  reduced  times  T = .01, 

.1,  1,  10,  as  labeled,  at  reduced  permeability  P = 0.5. 


17 


(x,T) 


Fig.  2(d).  Distribution  of  randomly  diffusing  particles  from  an 
instantaneous  source  at  the  center  of  the  middle 
layer  of  a 19-layer  system  at  reduced  times  T = .01, 
.1,  1,  500,  and  5000,  as  labeled,  at  reduced 
permeability  P = 0.001. 
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The  form  has  the  appearance  of  a Gaussian  or  of  segments  thereof. 
For  the  case  of  arbitrary  P and  very  short  T,  where  the  distribution 
function  is  small  at  the  edges  of  the  center  region,  or  for  the  case 


of  arbitrary  T but  infinite  P (no  barriers)  the  numerical  results 
quantitatively  fit  the  function 

p = 1/(4  nT)  ^ exp[-x  /4T] 


(10) 


In  the  case  of  zero  permeability  the  results  for  any  {x^,  T)  agree 

accurately  with  results  derived  earlier  from  a different  starting 

8 9 

formula.  ’ It  is  of  interest  to  note  that  for  this  limiting  case 
p may  also  be  accurately  constructed  by  repeated  total  reflection 
of  Eq.(10)  at  the  barriers. 

For  finite  permeability  and  diffusion  time  p is  not  Gaussian, 
nor  can  it  be  constructed  by  any  simple  scheme  of  partial  reflections 
of  a Gauss  function  at  the  barriers.  This  can  be  shown  from  the 
boundary  conditions,  Eqs.  (6)  and  (7). 

Finally,  it  was  found  that  the  numerical  values  at  the  midpoints 
of  the  regions  for  the  case  P = 0.001  do  not  follow  either  a Gaussian 
or  a binomial  distribution,  but  seem  to  approach  them  as  T inc»  xases. 

Diffusion  Coefficients 

Here  we  present  diffusion  coefficients  calculated  for  a wide 
range  of  barrier  permeabilities,  of  diffusion  times,  and  of  the 
gradient  strength  parameter  o Ycpsa.  With  the  19-region  system 
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considered,  it  was  possible  to  obtain  diffusion  coefficients  well 
into  the  range  at  long  diffusion  times  where  Oit)^  D(°°) , without 
violating  the  condition  that  the  effects  of  the  outer  boundaries  be 
smal 1 . 

In  Figure  4a  we  present  calculated  echo  heights  R(t)  versus 
reduced  diffusion  time  for  several  values  of  the  reduced  barrier 
permeability.  This  is  the  form  in  which  experimental  spin-echo 
diffusion  results  are  typically  presented.  For  low  barrier  perme- 
abilities ( P<1 ) it  is  the  most  useful  type  of  plot  for  deciding  whether 
the  experimental  diffusion  time  has  been  extended  to  long  enough 
values  so  that  D(™)  can  be  approximated. 

The  curves  all  approach  the  same  limiting  slope  as  T-O,  and  they 
appear  to  approach  various  smaller  limiting  slopes  as  T-**>,  though  we 
will  see  later  that  this  latter  is  true  only  for  o<l.  The  chord  to 
the  origin  from  any  point  is  proportional  to  the  apparent  diffusion 
coefficient  at  this  time  as  defined  by  Eq.(3).  For  a convenient 
comparison  these  diffusion  coefficients  are  plotted  on  the  same  time 
scale  in  Figure  4b.  When  = 0,  then  R becomes  independent  of  D(o) 
at  long  T,  and  is  determined  only  by  the  barrier  spacing,  according 
to8,9  R(°o)=  2(l-coso)/o  . 

In  Figure  5 the  diffusion  coefficients  are  plotted  versus  the 
logarithm  of  the  diffusion  time  in  order  to  cover  a larger  range  of 
the  latter.  A small  value  of  o is  chosen.  For  TO,  0(t)/D(o)  appears 


Fig.  4.  (a)  Semi-log  plot  of  echo  height,  R,  versus  reduced  diffusion 

time  for  o = 2.5,  at  the  reduced  permeabilities  indicated,  (b) 
Relative  diffusion  coefficients  on  the  same  time  scale  as  for 
the  parameters  of  (a). 


to  approach  unity,  as  expected.  At  long  diffusion  times  it  is  verified 
numerically  that  the  assymptotes  satisfy 

D(o)/D(»)  = 1 + 1/P,  for  o-*0.  (11) 

This  is  a relation  readily  derived  for  the  case  of  steady  state 
diffusion  from  one  side  to  the  other  in  a system  such  as  ours,  and 
where  D(°°)  would  represent  an  overall  diffusion  coefficient,  and 
D(o)  an  intrabarrier  diffusion  coefficient. 

A close  examination  of  the  two  curves  for  finite  permeability, 

P = 0.5,  2,  shows  that  they  are  not  even  approximately  linear 
combinations  of  the  curves  for  P=0,°°  . Therefore,  the  diffusion 
coefficients  for  a system  of  permeable  barriers  cannot  be  calculated 
by  a simple  superposition  of  results  for  impermeable  barriers  with 
results  for  no  barriers,  as  in  Eq.(l).  It  appears  that  the  drop 
in  diffusion  coefficient  takes  place  at  progressively  earlier  times 
as  the  permeability  increases.  Note  that  the  curves  for  P=2,<»  break 
away  at  the  longest  time,  where  the  effect  of  outer  boundaries  has 
become  important.  To  obtain  useful  results  in  this  region  we  would 
have  to  consider  a larger  number  of  layers. 

In  Figure  6 (solid  lines)  we  find  that  the  NMR  diffusion  coefficient 
as  defined  by  Eqs.(3)  and  (4)  is  not  solely  determined  by  the  system 
parameters  (D(o),  a,  £)  and  the  diffusion  time,  but  also  by  the 
strength  of  the  field  gradient  used  to  make  the  measurement.  Although 
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the  NMR  D(t)  becomes  independent  of  o and  also  agrees  with  the 
Einstein  D(t)  as  defined  by  Eq.(4)  when  0 approaches  zero,  D(tJ  for 
the  NMR  case  becomes  strongly  o-dependent  where  o>l.  This  is  true  for 


all  permeabilities  except  in  the  limit  P-*»  . The  variation  of  D(t) 
with  0 is  generally  monotonic  except  for  a few  small  regions  of  (0, 

T,  P)  space  where  either  R or  D increases  with  increasing  0. 

It  is  unfortunate  for  evaluating  experimental  results  that  the 
strong  o-dependence  occurs  in  a range  of  0 which  is  expected  for 
typical  experimental  systems,  such  as  those  of  Table  I.  The  fact 
that  the  values  of  D at  larger  0 go  below  the  expected  limits  at 
long  time  as  calculated  by  Eq .(11)  is  even  more  troublesome.  It 
is  physically  unreasonable  either  that  D should  increase  as  T is 
increased  beyond  the  range  presented  in  the  figures  or  that  D(«) 
should  be  less  than  as  calculated  from  Eq .(11). 

It  may  be  asked  whether  this  theoretical  result  is  an  artifact 
of  the  particular  geometry  chosen  or  of  its  ideal  regularity.  A 
satisfactory  answer  cannot  be  given  here,  and  awaits  calculations 
with  systems  of  other  geometries.  It  is  the  author's  belief  that 
it  is  the  results  for  small  0 which  will  more  likely  be  approximated 
when  calculations  with  other  geometries  are  performed,  and  that  these 
results  should  be  more  trustworthy  for  application  to  experimental 
systems  in  the  meantime. 
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TABLE  I 


Approximate  range  of  parameters  in  typical  pulsed  gradient  experiments 
on  colloidal  systems.  State-of-the-art  equipment  assumed  to  include 
capability  for  stimulated  echo  and  alternating  gradient  pulse  sequences, 
for  time  averaging  (used  at  long  t),  and  for  gradients  to  500  G/cm. 


Example 

Single  Cell 

Muscle  Cell 

Oil  Droplet 

System  parameters 

a , pm 

3 

50 

10 

D(o),  cm2/sec  xlO5 

0.5 

1.5 

0.01 

D(°°)/D(o) 

0.4 

0.6 

0 

T2,  msec 

40 

50 

50 

Tj,  msec 

500 

1000 

50 

Experimental  parameters 

t,  msec 

0.5  to  500 

0.5  to  2000 

1 to  60 

T 

0.03  to  30 

.0003  to  1.2 

.0001  to  .0006 

0 

1 to  0.2 

8 to  0.8 

17  to  10 

In  any  case  the  correctness  of  the  present  calculations  for  the 


particular  geometry  chosen  seems  assured  by  the  agreement  of  the 
results  for  £ and  D with  other  formulas  in  all  of  the  limiting  cases 
tried,  except  at  large  o and  T.  (See  Appendix  B for  additional  tests.) 

The  dashed  curve  of  Figure  6 corresponds  to  D(t)/D(o)  which  one 

2 

would  calculate  by  plotting  InR  vs  £ at  various  T,  and  estimating 
an  average  slope  over  the  first  decade  of  decrease  in  R to  calculate 
D.  Since  some  arbitrariness  was  involved  in  deciding  what  is  the  average 
slope  of  a curved  line,  this  dashed  line  is  only  approximate.  It 
is  presented  because  it  resembles  the  procedure  one  would  follow 
in  treating  real  data.  We  note  that  it  does  not  belong  to  the  family 
of  curves  represented  by  the  solid  lines.  This  is  because  in  actual 
experiments  we  tend  to  decrease  o as  we  increase  T,  so  as  to  prevent 
too  much  signal  attenuation,  whereas  the  solid  lines  are  at  constant  0 . 

The  times  at  which  D(tJ  reaches  its  arithmetic  mean  value  are 
plotted  in  Figure  7 for  each  of  several  0 . Times  for  the  arithmetic 
and  geometric  mean  values  at  small  0 are  presented  in  Table  2. 

Also  presented,  in  the  last  column  of  Table  2,  is  the  fraction  of 
particles  outside  the  center  layer  at  the  time  of  geometric  mean 
D(t).  The  significance  of  these  is  discussed  in  the  following  two 
sections. 
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TABLE  II 


Diffusion  times  and  particle  distributions  corresponding  to  the  mean 
of  the  initial  and  final  diffusion  coefficients 


At  Arithmetic 

Mean  D 

At  Geometric 

Mean  D- 

p 

D(-)/D(o) 

D(mean) 

T 

D(mean) 

T 

_/pdx 

5. 

.833 

.917 

,0081 

.910 

.0105 

i^o 

.063 

2. 

.667 

.833 

0257 

.815 

.035 

.088 

.5 

.333 

.667 

069 

.577 

.126 

.099 

.001 

.001 

.501 

113 

.032 

2.76 

.0056 

I 
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Physical  Interpretation 


It  was  initially  expected  that  the  time  of  mean  JD  would  not 
depend  strongly  on  P.  In  that  case  it  could  have  been  said  that 
the  observation  time  required  to  detect  the  presence  of  barriers 
depended  on  their  spacing,  but  not  on  their  permeability.  In  other 
words,  their  effect  would  be  proportional  to  the  fraction  of  particles 
which  had  contacted  them.  We  see  from  Figure  7 and  Table  2 that 
this  is  not  true. 

1 ? 

It  has  also  been  proposed  that  the  restrictive  effects  are 
manifested  in  proportion  to  the  fraction  which  have  been  reflected 
by  the  barriers.  In  that  case  the  time  of  mean  D(t)  should  increase 
with  increasing  P.  We  see  (Figs.  6 and  7)  that  just  the  opposite 
occurs. 

The  best  qualitative  description  found  is  that  the  effects  are 
manifested  in  proportion  to  the  number  of  particles  which  have 
penetrated  the  barriers.  This  is  illustrated  in  the  last  column  of 
Table  2 by  the  nearly  constant  value  of  the  fraction  of  particles 
outside  their  layer  of  origin  at  the  time  of  mean  D(t:)  for  all  but 
very  low  P.  Of  course,  such  a criterion  cannot  even  be  applied  to 
the  case  P = 0. 
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Evaluation  Of  Experimental  Results 


Thus  far  we  have  calculated  apparent  diffusion  coefficients, 
given  barrier  spacings  and  permeabilities,  and  intrabarrier  diffusion 
coefficients.  In  experimental  situations  we  desire  to  do  the  reverse— 
evaluate  the  barrier  parameters  and  intrabarrier  diffusion  coefficient 
from  a set  of  apparent  (measured)  diffusion  coefficients  over  a range 
of  diffusion  times. 

Although  the  actual  geometry  of  the  restrictive  barriers  may  be 
considerably  different  from  parallel  planes,  it  is  believed  that 
an  application  of  the  results  pertaining  to  small  o obtained  for  this 
case  should  give  a reasonably  good  approximation  to  the  correct  system 
pa rameters. 

Although  it  would  be  possible  to  fit  a set  of  data  to  Eq. (8) 
and  extract  a,  £ and  D(o)  by  a non-linear  regression  procedure,  this 
would  be  costly.  Instead,  it  is  worth  the  trouble  of  extending  the 
range  of  diffusion  times  as  far  as  possible,  so  that  £(o)  and  D(”°) 
can  be  estimated  directly.  The  product  ap  can  then  be  calculated 
from  Eq.(ll).  If  a is  not  known  from  other  information  it  may  be 
obtained  from  Figure  7,  employing  the  curve  for  small  0 . 


32 


SUMMARY 


A general  treatment  of  time-dependent  (transient)  diffusion 
coefficients  in  a system  containing  barriers  of  arbitrary  perme- 
ability has  been  given.  An  exact  expression  in  series  form  has  been 
obtained  for  the  distribution  function  at  arbitrary  time,  resulting 
from  one-dimensional  diffusion  originating  at  any  point,  and  occurring 
in  a homogeneous  medium  containing  equally  spaced,  planar  barriers 
of  arbitrary  but  equal  permeability.  This  function  has  been  used 
to  calculate  the  time  dependence  of  the  diffusion  coefficient  as  it 
would  be  measured  in  pulsed  gradient  NMR  experiments  employing  dif- 
fusion times  short  enough  so  that  few  of  the  diffusing  particles  had 
contacted  the  barriers,  to  times  long  enough  so  that  most  of  them 
had. 

The  time  dependent  diffusion  coefficient  was  also  calculated 
using  the  Einstein  formula  D(t)  = x7  /2t.  Numerical  results  of  the 
two  methods  agree  in  the  limit  where  the  applied  field  gradient  is 
zero.  However,  results  obtained  from  the  pulsed-gradient  formula 
are  strongly  dependent  on  the  gradient  when  y£<^a  exceeds  unity. 

Numerical  results  for  the  distribution  function  and  for  diffusion 
coefficients,  covering  a wide  range  of  reduced  diffusion  times, 
relative  permeabilities,  and  field  gradients  are  presented.  From 
the  results  it  is  shown  how  in  many  cases  the  reverse  process  of 
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evaluating  these  parameters  from  experimental  diffusion  data  by  hand 
calculation  is  possible. 


It  is  suggested  that  the  results  should  be  approximately  correct 
for  other  systems  of  repeated  barriers.  In  the  case  of  pulsed  gradient 
NMR  experiments  this  would  mean  colloidal  systems  such  as  biological 
tissues,  interconnected  pores  in  a rock,  or  the  continuous  phase  of 
an  emulsion. 
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APPENDIX  A 


NOTES  ON  THE  CALCULATIONS 
The  Eigenfunctions  of  p 

The  roots, 3 , of  the  matrix  of  boundary  equations  used  to 
determine  the  eigenfunctions  from  which  p is  constructed,  occur  in 
groups  of  m,  where  m is  the  number  of  layers  in  the  system.  Each 
group  begins  approximately  at  an  integral  multiple  of  n . The  spacing 
within  a group  increases  as  larger  values  of  P are  selected.  For 
P=-*>  we  have  3n=iWm. 

~ 2 

The  computer  time  required  is  approximately  proportional  to  m , 

because  the  number  of  3 required,  aid  the  time  required  to  evaluate 
the  sparse  matrix,  each  increase  linearly  with  m.  Most  of  the  time 
is  used  in  evaluating  the  s . Once  this  is  done  for  a give  P,  then 
p and  D(t)  (by  Eq.  (2)  or  (4))  may  be  calculated  quickly  for  a large 
number  of  locations  or  times,  respectively.  As  a typical  example 
for  a 19-layer  system,  the  time  required  to  evaluate  266  values 
of  3 (the  usual  number  of  terms)  to  double  precision  accuracy  (single 
precision  would  have  sufficed),  and  subsequently  evaluate  D by  Eqs. 

(2)  and  (4),  at  four  values  of  0 , each  at  70  values  of  T,  was  70 
sec.  on  a CDC  6600  computer. 

Listing  of  the  various  programs  are  available  on  request. 
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APPENDIX  B 


M 
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ADDITIONAL  CHECKS 

Several  additional  checks  on  the  correctness  of  the  results  are 
as  follows: 

It  was  found  that  R=1  for  0=0  for  all  T up  to  where  the  outer 
boundaries  become  important.  This  was  done  at  several  P. 

D{_t)/D(o)  was  calculated  directly  from  p by  a Simpson  rule  integration 
of  Eq.(2).  This  was  done  at  several  P and  several  T,  for  G>1 . 

The  results  were  identical  to  those  obtained  from  Eq.(8)  to  the  number 
of  significant  figures  printed  out. 

Equation  (8)  was  rederived  from  Eqs.  (2)  and  (5)  using  two  other 
coordinate  systems.  In  one  case  a single  origin  for  x and  was 
located  at  an  outer  boundary.  In  the  other  case  a single  origin 
was  located  at  the  center  of  the  center  layer.  Numerical  results 
were  in  agreement  with  those  for  the  coordinate  system  described 
earlier  for  all  values  calculated,  to  within  the  4-6  significant 
figures  printed  out. 

Distribution  functions  and  diffusion  coefficients  calculated 
for  P=0  agreed  to  the  number  of  figures  printed  out,  over  the  entire 

o q 

time  range,  with  the  results0’3  which  had  been  obtained  earlier  for 
this  special  case  using  a different  formula. 


A few  calculations  were  performed  with  a 41 -layer  system,  so 
as  to  extend  the  range  to  longer  T.  However,  D continued  to  decrease 
with  increasing  T for  all  0 used.  All  calculations  of  p or  D agreed 
among  systems  of  varying  number  of  layers,  provided  the  comparison 
was  at  times  short  enough  so  that  p was  negligible  in  the  outermost 
layers  of  the  smallest  system  considered. 
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